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ABSTRACT 

The observed correlation between global low cloud amount and the flux 
of high energy cosmic rays supports the idea that ionization plays a crucial 
role in tropospheric cloud formation. We explore this idea quantitatively with 
a simple model linking the concentration of cloud condensation nuclei to the 
varying ionization rate due to cosmic rays. Among the predictions of the model 
is a variation in global cloud optical thickness, or opacity, with cosmic-ray 
rate. Using the International Satellite Cloud Climatology Project database 
(1983-1999), we search for variations in the yearly mean visible cloud opacity 
and visible cloud amount due to cosmic rays. After separating out temporal 
variations in the data due to the Mt. Pinatubo eruption and El Nino/Southern 
Oscillation, we identify systematic variations in opacity and cloud amount due 
to cosmic rays. We find that the fractional amplitude of the opacity variations 
due to cosmic rays increases with cloud altitude, becoming approximately zero 
or negative (inverse correlation) for low clouds. Conversely, the fractional 
changes in visible cloud amount due to cosmic rays are only positively-correlated 
for low clouds and become negative or zero for the higher clouds. The opacity 
trends suggest behavior contrary to the current predictions of ion-mediated 
nucleation (IMN) models, but more accurate temporal modeling of the ISCCP 
data is needed before definitive conclusions can be drawn. 
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1. Introduction 



The primary source of energy for the Earth's atmosphere is the Sun, so it is reasonable 
to explore whether changes in the global climate result from solar variability. It was first 
suggested by the astronomer William Herschel ( |Herschel 1801| ) that variations in the solar 



irradiance caused by sunspots could lead to climatic changes on Earth, and he cited the 
variation of British wheat prices with sunspot number as evidence for this link. The 
occurrence of the "Little Ice Age" during the 1645-1715 Maunder sunspot minimum ([Eddy] 



1976|) , the correlation between the long-term solar cycle variations and tropical sea surface 
temperatures ( [Reid 1987|) , polar stratospheric temperatures ( [Labitske 1987| ), and the width 
of tree rings ( |Zhou and Butler 1998] ), along with many other studies also support a link 
between solar variations and the Earth's climate. 

A direct link between the Sun and these phenomena is tenuous, however, because the 
magnitude of the solar irradiance variation over the 11-year solar cycle is very small. Over 
the 1979-1990 solar cycle, for example, the variation in the irradiance was only ~ 0.1% 
( [Frohlich 2000 ), or ~ 0.3 W m" 2 globally -averaged at the top of the atmosphere. This is 



insufficient to power the sea surface temperature changes associated with the solar cycle by 
a factor of 3 — 5 ([Lean 1997| ) , and is significantly smaller than the globally-averaged forcings 
due to clouds (~ 28 W m -2 ; e.g. [Hartmann 1993|) , anthropogenic greenhouse gases (~ 2 W 
m -2 ; |Wigley and Raper 1992| ), and anthropogenic aerosols (~ 0.3 — 2.0 W m -2 ; ICharlsori 



et al. 1991% |Kiehl and Bricglcb 1993] ), suggesting that any direct atmospheric forcing from 



solar irradiance variations would be relatively unimportant. 

An indirect link between solar cycle variations and the Earth's climate appears more 
likely, especially given the discovery of a link between the flux of Galactic cosmic rays 
(GCRs) and global cloudiness (Svcnsmark and Friis-Christcnscn 1997) in the ISCCP cloud 



database ( |Rossow and Schift'er 19991) . The Sun modulates the GCR flux at the Earth 
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through the action of the solar wind, which scatters and attenuates the GCRs in times of 
heightened solar activity (solar maximum; e.g. |Jokipii 1971| ). Using 3.7 /iin infrared (IR) 



cloud amounts from the ISCCP database for the years 1983-1993, Marsh and Svensmark 
(2000) and Palle Bago and Butler (2000) showed that there is evidence of a positive 
GCR-cloud correlation only for low (< 3 km) clouds, and that the effect of the cosmic 
rays on global cloud amount appears to be greatest at the low to mid latitudes. The 
globally-averaged forcing due to the increase in low clouds associated with the solar cycle 
GCR variations is estimated ( Kirkby and Laaksonen 2000|) to be approximately —1.2 W 



m -2 , which is sufficient to power the sea surface temperature variations ( [Lean 1997|) . This 
is also comparable in magnitude (but opposite in sign) to the forcing due to anthropogenic 
CO2 emission over the last century (|Svensmark and Friis-Christensen 1997|) . Decreasing 
local cloud amounts correlated with short-term Forbush decreases in cosmic-ray rates were 
observed by Pudovkin and Veretenko (1995). 

The reality of the GCR-cloud connection has been questioned by a number of authors 
( [Kernthaler, Toumi, and Haigh 1999| ; LJorgensen and Hansen 2000) ; |Norris 2000|) . These 



objection can be distilled into three main points: 1) The GCR-cloud correlation should be 
seen prominently in high (cirrus) clouds at high latitudes where the cosmic-ray intensity is 
highest, 2) the increased cloudiness can be more plausibly attributed to other phenomena 
instead of GCRs, and 3) the correlation is an artifact of the ISCCP analysis. The first 
objection is addressed by the theory of ion-mediated nucleation (IMN: |Yu and Turco 2001] ; 



Yu 2002), in which the efficiency of the cosmic-ray interaction is limited at high altitudes 
by the lack of aerosol precursor vapors such as H2SO4 relative to the ion concentration. For 
the second objection, the temporal profile of the GCR-cloud correlation may be inconsistent 
with the profiles of the dominant volcanic and El Nino/Southern Oscillation (ENSO) events 
during the same time period ([Kirkby and Laaksonen 2000] ), although no quantitative study 
of the various temporal signatures in the data has been undertaken. Finally, the ISCCP 
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artifacts pointed out by Norris (2000) are troubling, but it is not clear that they are of 
sufficient magnitude to produce the observed GCR-cloud correlation, and it doesn't explain 
why the correlation exists only for low clouds and not the other cloud types in the ISCCP 
database. 

The linkage between cosmic rays and cloud formation has been recently investigated 
by a number of authors ( |Yu 2002] ; [Yu and Turco 20"0l| ; |Tinsley 2000| and references 
therein). Here we apply a perturbative approach to quantify the effects of variations in 
the cosmic-ray rate on the optical thicknesses, or opacities, of clouds, and use the observed 
cloud opacity variations to constrain the microphysical models of ion-mediated ultrafine 
particle formation. The paper is organized as follows. In the next section we discuss 
how the effect of cosmic rays could alter the optical thickness and emissivity of clouds by 
affecting the nucleation of condensation nuclei (CN). The search for variations in cloud 
optical properties using the ISCCP database and their correlation with cosmic-ray flux 
variations are discussed in Section ^] A discussion of the results is given in Section [4], and 
finally we summarize our results in Section |5] 



2. Effects of GCRs on Cloud Properties 



2.1. Nucleation 



Cosmic rays form water droplets in the supersaturated air of a classical cloud chamber 
( Wilson 19011) , and it seems plausible that they could also play a significant role in natural 
cloud formation. Yu and Turco (2000, 2001) and Yu (2002) have investigated the formation 
of ultrafine CN from charged molecular clusters formed from cosmic-ray ionization, and 
they find that the charged clusters grow more rapidly and are more stable than their 
neutral counterparts up to a size of ~ 10 nm. Although the subsequent growth of the 
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cosmic-ray formed ultrafine CN to viable ~ 100 nm cloud condensation nuclei (CCN) has 
not been explored, the concentration of CCN should also reflect the CN concentration, as 
well as the direct influence of cosmic rays, if the cosmic-ray ionization rate does not affect 
other important nucleation efficiency parameters such as condensible vapor concentration, 
temperature, and pressure. We will make this assumption here although it may not be 
strictly true with respect to the condensible vapor concentration (see e.g. [l'urco, Yu, and 
|Zhao 2000| ; |Yu 2002| ). 



Althought the formation of CCN and ultimately cloud droplets is a function of many 
variable factors such as temperature, pressure, vapor concentration, and relative humidity, 
we can quantify the effects of small variations in the ionization rate (primarily due to cosmic 
rays above ocean and at altitudes > 1 km above land; e.g. |Reiter 1992j ) on the number of 
CCN through a perturbation approach, i.e. 



O 7V J 

N CC K(q + Aq,V)^NcMq,V) + Aq ' ( ' ( " N 



(1) 

V 



dq 

where iVccN is the concentration of CCN, q is the ionization rate, V refers to the 
set of parameters other than the ionization rate affecting iVccN, and the partial 
derivative is evaluated for fixed V (hereafter this will not be written explicitly). Along 
with the assumption discussed previously, this approach assumes that the quantity 
Aq\dN cc ^/dq\ « Nccn(q,V), which is probably true for solar cycle variations, where q 
typically varies by < 30%, but may not be true during periods of large scale changes in the 
geomagnetic field (e.g. [lric et al. 1992] ). 



To quantify the effect of varying CCN concentrations on cloud optical thicknesses, we 
envision the two idealized scenarios depicted in Figure 1. In both cloud formation scenarios, 
changes in the ionizing cosmic-ray flux cause changes in the number of cloud condensation 
nuclei through the process of ion-mediated nucleation on the formation of ultrafine CN 
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in accordance with the assumptions mentioned above§. In the first case we assume 
that the nucleation of cloud droplets is limited by the available amount of water in the 
supersaturated air, so that the liquid water content (LWC), or density of water in droplets, 
is constant. Therefore the amount of water per droplet and hence the effective radii of cloud 
droplets will change with the cosmic-ray ionization rate. This is analogous to the "Twomey 
Effect" of enhanced aerosol pollution on droplet size distributions and the albedo of clouds 
( |l'womey 1977| ; |liosenfeld 2000[ ), and would primarily occur in environments where the 



amount of water in the air (and not the number of CCN) is the limiting factor. Thus, 
using ([J), we would expect that the effective radius R e ff of the cloud droplet distribution 
resulting from a small change in the cosmic-ray ionization rate Aq in any particular volume 
of air will be 

p [ N CCN (q,V) 1 1/3 / Aq <9iV CCN r 1/3 

6// " [N CCN (q + Aq,V)\ K "~[ 1+ N C ^ dq ) ^ 

where -Re// is the effective radius of the unperturbed droplet distribution, which we will 
associate with the solar maximum period of the solar cycle. 

In the second case in Figure 1, we assume that the change in CCN concentration 
resulting from change in cosmic-ray ionization causes a proportionate change in the amount 
of water extracted from the supersaturated air, with the effective radius of the cloud droplet 
distribution remaining constant. This is the case where the formation of the cloud is limited 
by the local availability of CCN and not condensible water. This effect has been seen in the 
marine boundary layer in ship track clouds ( |Conover 1966| ), which have higher reflectivities 
( |Coakley, Bernstein, and Durkee 1987] ) and liquid water contents ( [Radke, Coakley, and King] 



1989Q due to the formation of additional ultrafine CN from ship exhaust. The perturbed 



x We have assumed dN CCN /dq > in Figure 1, which need not be valid for all q. 
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liquid water content of a cloud in any particular volume of air will then be given by 

LWC « (l + LWCo, (3) 

V N CC N oq J 

where LWCo is the unperturbed cloud liquid water content associated with solar maximum 
as before. These two scenarios probably represent extremes of the direct cosmic-ray 
ionization effect on the clouds. As in the ship track clouds, the effect of the GCRs will 
probably be a combination of both LWC changes and R e ff changes, with the magnitude of 
the effect being bounded by the changes given in (g) and 

2.2. Radiative Properties 

Changes in the cloud liquid water content and droplet effective radius, associated with 
changes in the ionization rate due to cosmic rays, will result in changes in cloud opacities. 
The optical thickness r of a uniform cloud layer of thickness Az is given by ( |van den Hulst 



r = Az Q ext nfrWr dr, (4) 
Jo 

where n(r)dr is the concentration of cloud droplets with radii between r and r + dr, Q ext is 
the Mie extinction efficiency, and it is commonly assumed that 

Jo°° Qextn(r)r 2 dr _ 



/ °° n(r)r 2 dr 

which is a good approximation when 2nr/\ » 1, where A is the wavelength ( |Stephen"s 

mm. 



(5) 



The effective radius of the cloud droplet distribution is given by 



R = Jo 00 n(r)r 3 dr 
e ^ Jo°° n(r)r 2 dr ' 
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and the cloud liquid water content is given by 

4 r x 



4 rco 

LWC = -Tip / n(r)r 3 dr, (7) 
3 Jo 

where p is the density of liquid water. Combining these equations, we see that 

3 LWC Az 

Thus from (§) we would expect that an increase (decrease) in the mean R e ff and a decrease 
(increase) in the mean LWC, resulting from ionization variations due to cosmic rays, would 
result in a decrease (increase) the mean opacity of clouds. 

The change in cloud opacity with cosmic-ray rate can be quantified using the 
perturbation assumptions discussed in Section p.lj and equations (0), (|3|), and @. The 
fractional change in cloud opacity is then given by 

5r ^ Aq dN CCN 

r /iV CCN dq ' 1 ; 

where / = 1(3) for CCN (water) limited cloud formation, and the fractional change in the 
perturbed opacity r (relative to the unperturbed opacity r ) is defined by 5t/t = {t — tq)/tq. 
As mentioned previously, this derivation assumes that the right hand side of (§) is much 
less than one, which may not be the case for large changes in q and iVccN- As before we 
will assume that the unperturbed (perturbed) values of q and iV CCN refer to the values at 
solar maximum (minimum). 

At visible wavelengths from space, the primary consequence of the change in cloud 
opacity associated with cosmic rays will be an increase in cloud reflectivity, or albedo. To 
investigate this, we use the radiative transfer code SBDART QRicchiazzi, Yang, and Gautier 



to calculate the top of the atmosphere broadband (0.25 — 4.00 pm) upward flux for 



three uniform low cloud models: 1) a 1 km thick cloud layer extending to a height of 2 km, 
2) a 2 km thick cloud extending to a height of 3 km, and 3) a 0.5 km cloud layer extending 
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to 1.5 km. These simulations were done with a tropical atmosphere profile ( McClatchey et 



al. 1972|) and an ocean surface albedo. The fractional increases in albedo, resulting from 



a 10% increase in the number of cloud droplets due to cosmic ray ionization variations, is 
shown in Figure 2 for the 1 km thick cloud case, for a wide range of LWC and R e ff m the 
variable LWC case (top panel) and the variable R e ff case (bottom panel). In both cases 
the contours of changing albedo approximately parallel the change in optical thickness 
calculated assuming Q ext = 2.0 

Figure 3 shows the fractional change in albedo directly as a function of opacity for 
all three cloud models. This figure clearly shows that the change in albedo is largest for 
clouds with opacities r between 1 and 10, but is roughly independent of cloud geometrical 
thickness. Figures 2 and 3 indicate that the change in cloud optical thickness can be used to 
quantify the effects of the cosmic rays on cloud optical properties. Although the fractional 
change in albedo due to the cosmic rays is only ~ 2 — 5% for a 10% variation in the number 
of cloud droplets, this can produce a significant forcing per cloud of ~ 7 — 16 W m -2 at the 
top of the atmosphere for a solar zenith angle of 40°. The modulation of cloud opacity due 
to cosmic rays could therefore produce a similar modulation of the Earth's energy budget 
over the 11 year solar cycle, although the exact amount of forcing due to cosmic rays will 
depend sensitively on cloud amount variations, cloud opacity variations, and the efficiency 
at which changes in the cosmic ray rate are reflected in the number of cloud condensation 
nuclei. 

Because of the relationship between cloud opacity and emissivity, the cosmic rays 
should also produce an observable effect on cloud emission at infrared (IR) wavelengths. 
The effective IR emissivity e can be parameterized by a relation of the form ( |Stephen"s 



I978D 



e = 1 - exp(-a LWC Az), (10) 
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where a® is the mass absorption coefficient. Empirical fits to IR emission from water clouds 
yield a = 0.130 QStephens 1978|) . The exponent in ([U]) is proportional to the cloud optical 



thickness for a given droplet effective radius, so the infrared emissivity increases with cloud 
opacity, with the change being most noticeable for optically thin clouds. Therefore one 
would expect a change in IR emission, along with the primary effect of changes in visible 
albedo, from clouds at solar minimum relative to clouds at solar maximum if the cosmic 
rays change the cloud liquid water contents. Interestingly, a correlation between cosmic ray 
rate and cloud top temperature for low clouds has been reported by Marsh and Svensmark 
(2000), supporting this hypothesis. 

3. Cloud Opacity Variations 

3.1. ISCCP Data 



To search for systematic temporal changes in synoptic scale cloud optical properties, we 
used the International Cloud Climatology Project (ISCCP) monthly gridded cloud products 
( "D2" ) datasets, a compilation of cloud properties derived from satellite observations during 
the period 1983-1999 (|Rossow and Schiffer 1999| ). The ISCCP D2 data used here consists 
of mean daytime cloud amount fractions and visible optical depths, as a function of time, 
for 6596 "boxes" with equal area covering the entire surface of the Earth. For a given time, 
the cloud amount fraction in each box is defined as the number of cloudy satellite image 
pixels, as determined by a cloud detection algorithm, divided by the total number of pixels 
in the box. The cloud optical thicknesses are derived from the visible satellite cloud albedos 
by using a radiative transfer model and assuming spherical droplets with droplet sizes 
characterized by a gamma distribution with variance 0.15 and R e ff = 10 /im. ISCCP cloud 
top temperatures are simultaneously determined from the 3.7/xm IR radiances, allowing 
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for determination of cloud altitude and pressure, and the low, mid-level, and high clouds 
are defined as having cloud top pressures P > 680 mb, 440 < P < 680 mb, and P < 440 
mb, respectively. Because we require the simultaneous visible and infrared radiances to 
determine the opacity and cloud height for our analysis, we only use the ISCCP daytime 
data. This is a different dataset than the diurnal 1983-1993 IR data used for the cloud 
amount analyses of Marsh and Svensmark (2000) and Palle Bago and Butler (2000). 

Detailed information on the distribution of cloud optical thicknesses is not preserved in 
the the ISCCP D2 database, and instead the mean optical thickness fj is recorded for three 
broad opacity bands v. 0.0 — 3.6, 3.6 — 23.0, and 23.0 — 379.0. Thus a detailed analysis of 
the change in r over the solar cycle is not possible using the D2 data, but a value of the 
weighted mean cloud optical thickness f can be calculated using 

T= V » J , (11) 

2_i=i 

where the At are the total mean cloud amount fractions within each of the broad ISCCP 
optical thickness bins mentioned above. We calculated f separately for the three cloud 
altitude levels and for two latitude bands with \<j>\ < 40.0° (low latitude) and |0| > 40° (high 
latitude). The error associated with each fj was estimated by calculating the standard 
deviation of each ISCCP data point, from the scatter about the mean, and scaling by the 
square root of the number of data points. 

The mean optical thicknesses f as a function of time for the low latitude clouds are 
shown in Figure 4, and the corresponding result for global high latitude clouds is shown in 
Figure 5. Shaded is the two year period in which the effects of the Mt. Pinatubo eruption 
appear to be most significant. Also shown for comparison are the mean counting rates from 
the Climax, Colorado neutron monitor run by the University of Chicago (obtained from 
http://ulysses.uchicago.edu/NeutronMonitor/neutron_mon.html), which is a good measure 
of the local cosmic-ray ionization rate. In the low latitude case, the abrupt and large 
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decrease in f during 1991-1993 is due to the eruption of Mt. Pinatubo, and the subsequent 
plot scaling obscures smaller scale opacity variations. For comparison, we also plot the total 
mean cloud amount fractions A = Ya=i ^% for the same two latitude bands in Figures 6 and 
7. These plots show evidence for increases in mean cloud amount due to Mt. Pinatubo, as 
well as the smaller-scale temporal variations. 



3.2. Extracting The Cloud Variations Due to Cosmic- Rays 



To search for subtle variations in the ISCCP cloud opacities and amounts due to cosmic 
rays only, it is first necessary to eliminate the opacity variations in the data due to the Mt. 
Pinatubo volcanic eruption in June- September 1991 and strong ENSO events during the 
period of the ISCCP data. To separate out the various temporal signatures in the ISCCP 
data, we use a linear temporal model of the form 

3 

F(t) = J2hX k (t), (12) 

k=0 

where F(t) is the mean ISCCP quantity of interest for the year t, which for our purpose 
is either the visible cloud opacity f or the mean cloud amount /fractional area A. The 
model consists of four temporal basis vectors X^, which are functions of time, each scaled 
by a linear coefficient b^. For our temporal model we choose basis vectors corresponding 
to constant level of the given quantity (k = 0) and variations due to ENSO events (e.g. 
|Kuang, Jiang, and Yung 1998|) , the Mt. Pinatubo eruption of 1991, and cosmic rays 
(k = 1-3, respectively). Given the functional form of the basis vectors, the best-fit values 
of the linear coefficients can be determined through least squares minimization, and the 
fractional change in the time-varying ISCCP quantity over the data stretch is then given 
by 5F/F = bk/bo, where k = 1-3. This model assumes a linear correlation between the 
quantity of interest and the basis vectors and assumes no time delays; more complicated 
models are possible but will not be considered here. 
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The normalized basis vectors used in the temporal analysis of the ISCCP cloud 
data are shown in Figure 8. All of the vectors are scaled to values between zero 
and one. For the ENSO term X\ we use the scaled yearly-averaged Southern 
Oscillation Index (SOI) from the Australian Bureau of Meteorology (obtained from 
http: / / www.bom.gov.au/climate / current /soihtml .shtml|) . The SOI is a measure of the size 



of fluctuations in the sea level pressure difference between Tahiti and Darwin, Australia, 
and small values of the scaled SOI denote El Nino conditions and large values La Nina 
- both of which affect global weather ( [Rasmusson and Carpenter 1982 ). To parameterize 



the effect of the Mt. Pinatubo eruptions of 1991, we adopt a simple step function for X2, 
with identical non-zero intensities only for years 1991 and 1992. For the final term in the 
temporal model, X 3 , we use the scaled cosmic-ray rate from the Climax neutron monitor. 
Neutron monitor rates are directly proportional to the ionization rates due to cosmic rays 
because the neutrons are produced by the same cosmic ray cascade particles that produce 
the ionization, and the neutrons subsequently diffuse through less than 100 m of air before 
they are thermalized and captured by N to form 14 C (e.g. Lingenfelter, 1963). Neutron 
counters are thus unsusceptible to background ionizations due to terrestrial radiation from 
radioactive decays, which dominate the ionization signal from Galactic cosmic rays only 
below ~ 1 km in the atmosphere QReiter 1992| ). 



The results of the temporal fitting of both the ISCCP visible cloud opacities and 
amounts are shown in Table 1. Formally most of the fits are not good, with reduced 
chi-squares ranging from ~ 0.7 — 7.8 for twelve degrees of freedom. There are a number 
of possible factors that could be contributing to this. For example the error bars on the 
data may have been underestimated, leading to an artificially large values of chi-squared. 
Another possibility is that our fitting model is missing other significant temporal drivers, 
or perhaps a non-linear model or different basis vectors may be required to fit the data. We 
tried to fit the ISCCP data with linear models composed of different combinations of our 
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four basis vectors, and models with the cosmic-ray term provided a better fit to the data 
in general. Nevertheless it is possible that un-modeled phenomena mimic the temporal 
signature of cosmic rays in the data; more robust calculations of ISCCP error bars, inclusion 
of more ISCCP data, and exploration of more complicated temporal models in future work 
will help resolve this issue. 

The fractional variation in visible opacity 5t/t associated with the cosmic rays ranges 
from ~ +10% for high clouds to —7% for low clouds. For the mean visible cloud amounts 
the variation due to cosmic rays is just the opposite - becoming greater in magnitude as the 
cloud height decreases - qualitatively consistent with the positive correlation seen in the 
ISCCP IR data between cosmic ray rate and low clouds QSvensmark and Friis-Christcnscn 



1997| ; Marsh and Svensmark 2000; Palle Bago and Butler 2000). Therefore the high clouds 
appear to become thicker but smaller in response to increasing cosmic ray flux, while for 
the low clouds the response is the opposite. 

4. Discussion 



The observed variation of cloud optical thicknesses with cosmic-ray rate can be used 
to constrain microphysical models of the cloud condensation nuclei concentration Nqcn 
using (|9|). Of crucial importance is the partial derivative <9iV CCN /<9g, which determines the 
sign of the change in opacity with cosmic-ray rate. Recently Yu (2002) calculated iVccN 
as a function of altitude and ionization rate using an ion-mediated nucleation code. Given 
this model and the vertical profiles of sulfuric acid vapor concentration, ionization rate, 
temperature, relative humidity, pressure, and surface area of pre-existing particles assumed 
therein (Yu 2002), the value of dN CCN /dq peaks at values of q pea k = 12, 8, and 4 ion 
pairs cm -3 for low, mid-level, and high clouds, respectively, such that dN cc ^/dq > for 
q < q P eak and dNccTss/dq < for q > q pea k- Using the cosmic-ray ionization rates found by 
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Neher (1961,1967) interpolated to geomagnetic latitude 40°, we find typical ionization rates 
of, respectively, q ~ 3, 8, and 23 ion pairs cm -3 for the low, mid-level, and high ISCCP 
clouds. Therefore from @ we would expect a positive or zero correlation between opacity 
and cosmic-ray rate only for low clouds, and negative correlations for higher clouds for this 
model. We observe just the opposite, but the precision of the temporal model fits to the 
ISCCP data is not sufficient for us to rule out the Yu (2002) model based on the data. 

All three of the time-varying parameters in our temporal model show inverse 
correlations between mean visible cloud opacity and amount, suggesting a common origin 
for this behavior. These inverse correlations are illustrated in Figure 9. These are 
probably not artifacts of the averaging process because the quantities in Figure 9 have been 
normalized by their constant model terms in their temporal fits. One possible explanation 
for the inverse opacity-amount correlation is via a feedback mechanism. For the case of 
positive opacity variations, an increase in mean cloud opacity and albedo would result in 
increased energy loss to space and eventually less surface heating and subsequent water 
evaporation. Hence clouds would tend to be smaller and have smaller areas than they 
would otherwise. Conversely, for negative opacity variations clouds would tend to be 
larger. Global climate simulations ( |Chen and Kamaswamy 1996 ) indicate that global cloud 



albedo-increasing perturbations - similar to the changes induced by cosmic rays - decrease 
the global transport of moisture from the tropics, which then could conceivably produce 
fewer or smaller global clouds on average by this mechanism. Dynamical simulations of the 
response of global cloudiness to synoptic changes in the opacity are needed to investigate 
this. 



5. Summary 
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Here we consider a model in which Galactic cosmic rays alter the optical properties 
of clouds by changing the number of available cloud condensation nuclei. The main 
observational consequence of our model is a change in mean cloud opacity, with a secondary 
effect being a change in infrared emittance for optically thin clouds due to the relationship 
between cloud emissivity and opacity. We use the global ISCCP cloud database to search 
for variations in cloud properties due to cosmic rays, and after subtracting the background 
signals in the data due to Mt. Pinatubo and ENSOs, we find systematic variations in both 
opacity and cloud amount associated with changes in the cosmic-ray rate. The fractional 
variation in opacity attains a maximum positive value for high clouds and decreases with 
height, becoming negative or zero for low clouds. The fractional variation of the cloud 
amounts with cosmic-ray rate, however, show the opposite trend - increasing from a 
negative correlation at high altitudes to a positive correlation at low altitudes, which is 
consistent with the positive correlation between global low clouds clouds and cosmic-ray 
rate seen in the infrared ( |5vensmark and Friis-Christcnsen 1997 ; Marsh and Svensmark 



2000; Palle Bago and Butler 2000) 

Clearly more work is needed to model the opacity and cloud amount variations seen in 
the ISCCP data. Using our simple temporal model and perturbative approach, we have 
outlined a framework on which the variations in the data due to cosmic rays can be isolated 
and compared to model predictions. As the time span of the ISCCP data increases in 
length, more complicated models with additional components and nonlinear dependencies 
can be used, and the analysis can then be more robust. The ISCCP data requires the 
culling together and normalizing of many disparate satellite datasets flRossow and Schift'cr 



1999 ), and although this approach is necessary at the present time it is not ideal. One 
complement to the ISCCP global cloud data would be provided by the NASA deep space 
mission Triana, which would be able to retrieve cloud optical thicknesses simultaneously 
over the entire sunlit Earth from the LI Lagrangian point between the Earth and the Sun. 
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Continuous deep space observing of Earth's clouds would be ideal for detecting not only 
the solar cycle variations seen here but also the shorter duration but possibly more frequent 
variations in global cloud cover associated with Forbush decreases of Galactic cosmic rays 
and high energy solar proton events from the Sun. 

We thank the AVANTI article service of the Scripps Institution of Oceanography 
Library, and acknowledge the use of cosmic-ray data from the University of Chicago 
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Table 1. Temporal Fits to ISCCP Mean Cloud Opacity (f ) and Cloud Amount (A) 



Par. 


Lat. 


Alt. 


r 6 ° n 

[Const.] 


h 
[SOI] 


b 2 

[Pmatubo] 


[CR] 


■f 


< 40° 


Hif?h 


13 31 ± 23 


7.6 ± 1.7 


-32 2 ± 1 6 


9.4 ± 1.6 




< 40° 


Mid- level 


10 45 ± 32 


7.9 ± 2.9 


-28 3 ± 2 7 


5.8 ± 2.9 


j- 


< 40° 


Low 


3 78 ± 14 


-1.9 ± 3.1 


-4.2 ± 3.6 


-2.5 ± 3.4 


f 


> 40° 


High 


12.00 ±0.50 


16.6 ±4.0 


-6.9 ±3.6 


10.9 ±3.5 


f 


> 40° 


Mid-level 


6.81 ±0.23 


12.4 ±3.0 


-9.7 ±3.0 


1.8 ±3.2 


f 


> 40° 


Low 


5.01 ±0.14 


1.8 ±2.3 


-7.9 ±2.6 


-6.7 ±2.6 


A 


< 40° 


High 


13.84 ±0.52 


-4.8 ±0.32 


61.6 ±5.8 


-9.3 ±3.6 


A 


< 40° 


Mid-level 


5.71 ±0.16 


2.0 ±2.5 


58.1 ±4.2 


-0.1 ±2.7 


A 


< 40° 


Low 


45.5 ± 1.6 


6.9 ±3.3 


13.6 ±4.0 


0.8 ±3.5 


A 


> 40° 


High 


15.90 ±0.58 


-12.3 ±3.3 


7.9 ±4.0 


-9.9 ±3.3 


A 


> 40° 


Mid-level 


23.61 ±0.84 


-4.6 ±3.0 


8.3 ±3.8 


-3.1 ±3.4 


A 


> 40° 


Low 


38.7 ± 1.0 


7.1 ±2.4 


10.1 ±2.9 


6.5 ±2.6 



Except for the amplitude of the constant term 6 , the model amplitudes 
for the SOI, Mt. Pinatubo, and cosmic-ray terms (61-63, respectively) are 
normalized to 60 and expressed as a percentage. The mean cloud opacity f 
and total mean amount fraction A are from the ISCCP visible band data. 
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Fig. 1. — Cartoon illustrating two limiting scenarios for the effect of the Galactic cosmic 
rays (GCRs) on cloud optical properties, assuming that the varying ionizing cosmic-ray flux 
causes changes in the number of cloud condensation nuclei (CCN) through ion-mediated 
nucleation. In the first case we assume that the nucleation of cloud droplets is limited by 
the available amount of water in the supersaturated air. Therefore, as illustrated in the 
top panel, if an increase in the GCR ionization flux resulted in more cloud condensation 
nuclei (CCN) but no additional water condensation, the amount of water per droplet will 
be less and the effective radius R e ff of the droplet distribution will be smaller. Alternately, 
as illustrated in the bottom panel, if the formation of cloud droplets is limited by the local 
availability of CCN and not condensible water, R e ff can remained unchanged and additional 
CCN resulting from changes in cosmic-ray ionization would cause an increase in the amount 
of water extracted from the supersaturated air, so the cloud liquid water content would 
increase. The opposite trends hold for cases where the number of CCN is decreased by 
variations in the cosmic-ray flux. 

Fig. 2. — The fractional change in the albedo of a 1 km thick cloud expected from a 10% 
increase in the number of cloud droplets due to changes in the cosmic-ray flux, shown for 
the case of variable cloud water content LWC (top) and for variable droplet radius R e ff 
(bottom). The solid contours denote the change in albedo, and the dotted contours are for 
the optical thickness. 

Fig. 3. — The fractional change in the albedo, expected from a 10% increase in the number 
of cloud droplets from variations in the cosmic-ray rate, plotted as a function of cloud optical 
thickness for three different cloud geometrical thicknesses. The open symbols denote changes 
in cloud LWC and the filled symbols changes in R e ff- 
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Fig. 4. — The mean cloud 0.6 /zm optical thickness from the ISCCP database for all clouds 
in the low latitude band \(f>\ < 40°, with the cosmic-ray rate from the Climax neutron 
monitor. The high, mid-level, and low clouds refer to cloud top pressures of P < 440 mb, 
440 < P < 680 mb, and P > 680 mb, respectively, and the shaded interval refers to cloud 
data affected significantly by the eruption of Mt. Pinatubo in June 1991. The la error bars 
on the ISCCP data were calculated from the sample variance of the data. 

Fig. 5. — Same as Figure 4, but for all high latitude clouds with \(f>\ > 40°. 

Fig. 6. — The mean cloud 0.6 /im amount fractions from the the ISCCP database for all 
clouds in the low latitude band \(f>\ < 40°, with the cosmic-ray rate from the Climax neutron 
monitor. The la error bars on the ISCCP data were calculated from the sample variance of 
the data. 

Fig. 7. — Same as Figure 6, but for all high latitude clouds with \<j>\ > 40°. 

Fig. 8. — Basis vectors used in the temporal model of ISCCP visible opacity and amount 
fraction variations. The vectors X , X 1 , X 2 , and X 3 represent the constant level term and 
variations due to ENSO, the eruption of Mt. Pinatubo, and cosmic rays, respectively. 

Fig. 9.— Fractional change in visible cloud amount versus fractional change in visible 
opacity, from the fit of the temporal model to the ISCCP data. The points have been 
normalized by their respective constant term values in the temporal model, and the error bars 
have been omitted. In all cases there is an inverse correlation between the two parameters. 
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